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One of the standard generalized-gradient approximations (GGAs) in use in modern electronic- 
structure theory, PBE, and a recently proposed modification designed specifically for solids, PBEsoI, 
are identified as particular members of a family of functionals taking their parameters from difi'erent 
properties of homogeneous or inhomogeneous electron liquids. Three further members of this family 
are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules and 
solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, 
is not always the best performing member of the family, and that PBEsol, in spite of having been 
constructed specifically for solids, is not the best for solids. The performance of GGAs for finite 
systems is found to sensitively depend on the choice of constraints steaming from infinite systems. 
Guidelines both for users and for developers of density functionals emerge directly from this work. 
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Modern density-functional theory (DFT)ii^i^ owes its 
success and popularity largely to the availability of sim- 
ple and reliable density functionals 4. Among the most 
widely used such functionals are gradient-dependent ap- 
proximations, such as the B88 exchange functional^ or 
the PBE generalized-gradient approximation (GGA) for 
exchange and correlatioui^ Although many other func- 
tionals are available, PBE is today the de facto standard 
for gradient-dependent functionals in solid-state physics, 
and, together with B88, in quantum chemistry. These 
gradient-dependent functionals also form the basis for 
the development of more sophisticated functionals, of, 
e.g. the meta-GGA or hybrid type. 

Given the importance of PBE both for countless prac- 
tical applications of DFT, as well as for constructing 
more refined approximations, it is not surprising that 
over the years many variations of the basic PBE form 
have been develop ed J^i^ i ^ ° i ^ ^ Most of these are more em- 
pirical than the original construction, in the sense that 
they include parameters fitted to test sets of selected sys- 
tems and properties. None is uniformly better than the 
original PBE for all systems and properties. For more re- 
stricted classes of systems, however, it is not that hard to 
improve on PBE, as is illustrated by the recent proposal 
of PBEsolfi^ which was designed to improve on PBE 
specifically for solids. As this is still an extraordinary 
large and diverse class of systems, even an improvement 
'only' for solids is still a very major step forward, and 
consequently PBEsol is currently being implemented in 
many standard electronic-structure codes, and intensely 
scrutinizedii^ii^ 

In the present paper we point out that the step that 
led from PBE to PBEsol is, in fact, not unique, and al- 
lows several variations. We propose a family of function- 
als, which we call PBE(/3,/x), of which both the original 



PBE and the original PBEsol are particular members, 
and which includes at least three more alternatives. Each 
member of this family takes the value of the /3 and fi 
parameter from a different physical constraint, without 
statistical fitting to test sets of systems. Five PBE(/3,/Lt) 
functionals are systematically tested for atoms, molecules 
and solids, and, within each class, for systems with phys- 
ically and chemically different properties. 

To begin, we recall that the original PBE functional 
has three parameters. One, /3, appears in the correla- 
tion functional, and was originally obtained by requir- 
ing that for weakly inhomogeneous high-density systems 
the second-order gradient expansion of the correlation 
energy was recovered?!-- Another, /i, appears in the ex- 
change functional, and in the original construction was 
obtained from requiring that the combined exchange 
and correlation functional predicts the correct linear- 
response of bulk Jelliumi^ This latter requirement im- 
plies fi = 7r^/3/3. In PBE, this is used to fix /i, once f3 is 
obtained from the correlation energy gradient expansion. 
(One of the alternatives we propose here is to invert this 
relation, using it to determine /3 once /i is obtained from 
the gradient expansion of the exchange energy.) 

In PBEsol, (3 is determined instead by fitting to Jel- 
lium surface energies (JSEs) (a strategy proposed earlier 
by Armiento and Mattsson^^), whereas /i is chosen such 
as to recover the second-order gradient expansion of the 
exchange functional. Our original motivation for recon- 
sidering the choice of (3 and /i was that the JSEs used for 
obtaining /3 in PBEsol had themselves been calculated 
from the TPSS meta-GGA functionaljii whose correla- 
tion energy in turn contains as an ingredient the original 
PBE correlation. It seems slightly inconsistent to us to 
employ in the construction of a functional designed to im- 
prove on PBE, a fit to data obtained from a functional 
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TABLE I: Specification of tfie values of /3 and /i used in eacfi of tfie investigated variants of PBE(/3, fi). 

PBE(Gc,Jr) PBE(J.,G:.) PBE(J.,JO PBE(Gc,G'.) PBE(J.,G.) PBE(Gc,Js 

= PBE = PBEsol new new new new 

P 0.066725 0.046 0.046 0.066725 S^/tt^ 0.066725 

^ 71^/3/3 10/81 TT^/J/S 10/81 10/81 n.a. 



that itself contains PBE. In practice, however, this may 
be a rather purist objection, as the TPSS JSEs are quite 
close to those obtained by other methods f^^ii? 

A search for possibilities to avoid this small inconsis- 
tency, however, suggested a broader perspective on PBE 
and PBEsol: since the aim of PBEsol was to improve on 
PBE specifically for solids, and the gradient expansion is 
more relevant for the slowly varying densities of typical 
solids than it is for molecular densities, it seems promis- 
ing to take both parameters, /3 and fj,, from the gradient 
expansion, without making any use of properties of Jel- 
lium. An additional advantage of extracting both param- 
eters from the same source is that it enhances the chances 
of error compensation between the exchange and the cor- 
relation functional - an effect that is known to be behind 
the success of the local-density approximation, and that 
has only partially been preserved by common GGAs. 

Further exploration of the idea of a consistent (in the 
sense of coming from the same source) set of parameters 
suggests the alternative possibility to take both parame- 
ters, /3 and /u, from properties of Jellium, without making 
any use of gradient expansions. Specifically, we can de- 
termine f3 from JSEs and /i from the Jellium response 
function. A priori, this choice, too could be expected to 
be good for solids, or at least for metals, as Jellium is the 
paradigm of metallic behavior in extended systems. 

Table U presents a summary of the different variations 
of PBE and the corresponding values of /? and /i. Our 
notation for the entire family of functionals is PBE(/3, /j,), 
where f3 and are replaced by Gx (Gc) or Jg {Jr), de- 
pending on whether the parameter has been obtained 
from the gradient expansion for exchange (correlation), 
or from Jellium surface (response properties) . In this no- 
tation the original PBE becomes PBE(Gc, Jr), the orig- 
inal PBEsol becomes PBE( Js, Gx), and the two alterna- 
tives just described read PBE(G'c, Gx) and PBE{Js, Jr)- 

Table [J also includes a further mixed choice, in which 
(3 is determined from the Jellium response function and fi 
from the gradient expansion of the exchange energy, i.e. 
PBF,{Jr,Gx). This is the functional alluded to above, 
in connection with inversion of the relation /i = 7r^/3/3. 
Mathematically, still another possibility would be to de- 
termine P from the gradient expansion of the correla- 
tion energy and ^, e.g., from fitting to JSEs, resulting in 
PBE (Gc, Js), but the results from that fit are not avail- 
able, so we cannot test that option here. 

We note that PBE and PBEsol contain a third param- 
eter, K, which is determined from the Lieb-Oxford (LO) 
bound^S. Recently it has been suggested that the value 



TABLE II: Mean error (me), mean absolute error (mae) and 
mean absolute relative error (mare) (same labels are used in 
all Tables) for ground-state energies (GSE) Eq (Ry) and ion- 
ization potentials IP (Ry) for 26 light atoms^ (Z = 3 - 28); 
and relative percentage error for the H-atom GSE (RE^r). 







Gc , J r 


J S 1 Gx 


Js ^ Jr 


Gc-, Gx 


J r 1 G X 




me 


9.052 


7.790 


8.231 


7.594 


7.893 




mae 


9.354 


9.210 


9.232 


9.207 


9.211 




mare 


0.00451 


0.00627 


0.00560 


0.00662 


0.00609 


REh 


% 


0.00191 


2.26 


1.43 


2.65 


2.03 




me 


-0.245 


-0.247 


-0.244 


-0.254 


-0.243 


IP 


mae 


0.258 


0.260 


0.257 


0.266 


0.257 




mare 


0.368 


0.370 


0.367 


0.380 


0.366 



of this parameter can be readjusted to possible tighter 
forms of the LO bound. ^^'^^ In the present work we do 
not change k relative to the original PBE proposal, in 
order not to mix the question of a possible tightening of 
the LO bound with that of consistent choices of the /3 
and [I parameters. Future work should explore the con- 
sequences of combined changes of all three parameters. 
The proposal of Ref. [l^, which in addition to values of 
parameters such as k, /3 and /i also changes the form of 
the PBE enhancement factor, is a step in this direction. 

Restricting ourselves to exploring changes of only (i 
and /i, we report in Tables [iTl to IIVI the mean error, de- 
fined as me — jj X^iLi i^i ~ the mean absolute er- 
ror (mae — jj X^ili ~ ^^'^ mean absolute 
relative error {mare = ^f^i \{Ri ~ Ai) / Ri\) , where 
Ri is the reference (benchmark) value of the i'th system 
and Ai the corresponding approximate value, for TV = 26 
light atomsr^ = 14 molecules^! and A = 13 extended 
systems,— chosen according to two requirements: (i) reli- 
able benchmark data are available, and (ii) different types 
of crystal structures and of chemical bonds are repre- 
sented in the data set. In the following analysis, we focus 
mainly on the trends of the marc. 

For atoms, we have performed all-electron calculations 
using the mesh-based (basis-set free) atomic code that 
is part of the Siesta package.^ We find that the original 
PBE predicts best ground-state energies. The alternative 
functional PBE(Js, Jr), which takes both of its parame- 
ters from Jellium properties, comes second. This is rather 
unexpected, as JSEs (used to obtain p) and bulk Jellium 
response functions (used to obtain fi) are maximally dif- 
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TABLE III: The me, mae and mare for atomization ener- 
gies AE (eV) and interatomic distances d (A) for 14 small 
molecules.— 





PBE(/3,m) 


Gc 5 J r 


J S 7 Gx 


Jsi Jr 


Get Gx 


J r 5 Gx 




me 


0.0452 


0.566 


0.532 


0.458 


0.615 


AE 


mae 


0.299 


0.614 


0.584 


0.616 


0.647 




mare 


0.0573 


0.100 


0.105 


0.112 


0.103 




me 


-0.0143 


-0.0122 


-0.0115 


-0.0124 


-0.00959 


d 


mae 


0.0231 


0.0214 


0.0201 


0.0307 


0.0223 




mare 


0.0156 


0.0152 


0.0146 


0.0204 


0.0156 



TABLE IV: The me, mae, and mare for lattice constants a 
(A), bulk moduli B (GPa) and cohesive energies CE (eV) for 
13 solids.i^^ 





PBE(/3,m) 


G c , J r 


J S 5 Gx 


J Si Jr 


Gc, Gx 


J r 1 G X 




me 


-0.129 


-0.0722 


-0.0881 


-0.0692 


-0.0724 


a 


mae 


0.129 


0.0739 


0.0890 


0.0692 


0.0753 




marc 


0.0283 


0.0165 


0.0198 


0.0153 


0.0168 




me 


35.81 


25.65 


29.57 


24.06 


26.80 


B 


mae 


35.99 


25.84 


29.76 


24.21 


27.03 




mare 


0.244 


0.173 


0.201 


0.153 


0.184 




me 


-0.521 


0.00831 


-0.172 


0.0136 


-0.0136 


CE 


mae 


0.555 


0.282 


0.280 


0.288 


0.278 




mare 


0.117 


0.0604 


0.0659 


0.0644 


0.0598 



ferent from what one expects for atomic densities. Prob- 
ably, this indicates a substantial error cancellation be- 
tween exchange and correlation, made possible by taking 
both parameters from the same reference system. 

For atoms we also compared Koopman's theorem 
ionization potentials, i.e. the negative of the eigen- 
value of the highest occupied Kohn-Sham orbital. Here 
original PBE, PBE(J5, J^) and the novel mixed choice 
PBE( Jr, Gx) perform very similarly, whereas PBEsol and 
PBE(G'c, Ga;) are clearly worse. However, as ionization 
energies obtained from eigenvalues suffer from the self- 
interaction error and the resulting wrong asymptotics of 
the effective potentials, which is not addressed at all by 
changing the values of /3 and /i, a comparison of the per- 
formance of PBE(/3, [i) for these quantities is less conclu- 
sive than that for total energies. 

As a final test for atoms, we compared the error in the 
ground-state energy of the hydrogen atom. This error, 
arising in an A'^ = 1 electron system, is exclusively due 
to self-interaction. As the fifth line of Table |TT] shows, 
the original PBE functional performs by far best of all 
tested variants, distantly followed by PBE(J2,, J^)- As all 
constraints investigated here stem from infinite systems, 
this illustrates a dramatic sensitivity of the performance 
of approximate functionals near A^ = 1 to changing con- 
straints arising from N oo. 

For molecules and solids, our calculations were per- 



formed with the SIESTA code^^S employing norm- 
conserving TrouUier-Martins pseudopotentials and a 
strictly localized TZP and DZP basis set for molecules 
and solids, respectively. Pseudopotentials for use with 
any of the PBE(/3, /i) functionals were generated by em- 
ploying the same functional also in the atomic all-electron 
calculations. Convergence with respect to the basis set 
size was checked by repeating the molecular PBE calcu- 
lations, as well as selected atomic and solid-state calcula- 
tions, independently with a plane- wave code that is part 
of the CPMD packager^ where convergence is controlled 
by a single parameter, the energy cutoff. 

For molecules, original PBE=PBE(G'c, Jr) predicts 
best atomization energies, whereas the alternative 
PBE( Js,Jr) predicts best interatomic distances. For 
both quantities PBEsol=PBE( Jg, Gj,) comes second. 

For solids, the original PBE functional is the worst of 
all tested members of the PBE(/3, /i) family. Not surpris- 
ingly, PBEsol consistently outperforms PBE, but inter- 
estingly, it is itself outperformed by PBE(Gc, Gx) for lat- 
tice constants and bulk moduli, and by the novel mixed 
choice PBF,{Jr,Gx) for cohesive energies. Overall, the 
best performance is obtained from PBE(Gc, Gx)^^^ 

The present results suggest guidelines for the informed 
choice of density functionals and their construction: (i) 
The original PBE is not optimal for use in solid-state cal- 
culations and for molecular bond lengths, (ii) PBEsol, in 
spite of its name, is not the best PBE-type functional for 
bulk properties of solids, where PBE(Gc,Ga;) performs 
much better. Surfaces require separate analysis, (iii) A 
consistent choice of parameters (in the sense of stemming 
from the same type of source) tends to benefit the perfor- 
mance of the functional: PBE(Gc,Ga;) and PBE{Js, Jr) 
are more often among the best or second-best performers 
than any of the other variants. "^^ (iv) The good perfor- 
mance of PBE(Gc,Ga;) for solids suggests that restor- 
ing the gradient expansion completely (for exchange and 
correlation, and not just partially, for exchange only, as 
in PBEsol) is beneficial for extended systems, whereas 
for finite systems (atoms and molecules) gradient expan- 
sions turns out to be a less important ingredient, (v) 
A common ingredient of functionals that perform well 
for finite systems is the Jellium response function, but 
there seems to be no systematic trend for the utility of 
the Jellium surface energy, (vi) Substituting a constraint 
arising from infinite systems by another, also arising from 
infinite systems, can have a dramatic effect even at the 
opposite end of the spectrum, for A^ = 1. 

Users of PBE and related functionals should thus be 
aware of the fact that within the functional form of PBE 
there is not one single choice of constraints that bene- 
fits all system types. A related aspect is that the orig- 
inal PBE has been used as ingredient in the develop- 
ment of recent meta-GGA functionals, such as TPSSii 
and PKZBi^ These constructions should also be recon- 
sidered in view of the superior performance of some of 
the other PBE(/3, fi) variants, as more sophisticated high- 
level functionals should naturally be build using as ingre- 
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